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Abstract 

Some techniques for the use of bitwise operations are described 
in the article. As an example, an open problem of isomorphism-free 
generations of combinatorial objects is discussed. An equivalence 
relation on the set of square binary matrices having the same number 
of units in each row and each column is defined. Each binary matrix 
is represented using ordered n-tuples of natural numbers. It is shown 
how by using the bitwise operations can be implemented an algorithm 
that gets canonical representatives which are extremal elements of 
equivalence classes relative to a double order on the set of considered 
objects. 



1 Introduction 

The present study is thus especially useful for students educated to become 
programmers as well as for their lecturers. A meaningful example for the 
advantages of using bitwise operations for creating effective algorithms in 
programming is presented in this article. We will consider an open combi- 
natorial problem on binary matrices and its solution using the algorithm for 
some values of the integer parameters n and k. To implement the algorithm, 
we will use essentially bitwise operations. 

The use of bitwise operations is a powerful method used in C/C++ and 
Java programming languages. Unfortunately, in the widespread books on 
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this topic there is incomplete or no description for the work of the bitwise 
operations [21 [TJ El EH HSl CI]- The aim of this article is to correct this lapse 
to a certain extent and present a meaningful example of a programming 
task, where the use of bitwise operations is appropriate in order to facilitate 
the work and to increase the effectiveness of the respective algorithm. 
This work is an extension and complement to [S]. 

A binary (or boolean, or {0,l)-matrix) is a matrix whose all elements 
belong to the set B = {0, 1}. With Bn we will denote the set of all ri x n 
binary matrices. 

Some algorithms for isomorphism-free generations of combinatorial ob- 
jects are discussed in detail in [H]. In our work we will consider a problem 
of this type. Its formulation is as follows: A set of binary matrices C O Bn 
is given. In C is defined an equivalence relation. An algorithm which did 
not study every element of the set C, and which receives one representative 
of each equivalence class to be described. For this purpose, we will use 
significantly bitwise operations. 



2 Task formulation 

Let n and k be positive integers. Using the notation from [15], we will call 
A^-matrices all n x n binary matrices in each row and each column of which 
there are exactly k in number I's, and with A^ we will denote the set of 
these matrices. 

Let A,B e A^. We will say that A ~ 5, if A is obtained from B by 
moving some rows and/or columns. Obviously, the relation defined like that 
is an equivalence relation. We denote with 



li{n,k) 

the number of equivalence classes on the above defined relation. 

Problem 1 Find fi{n, k) for given integers n and k, 1 < k < n. 

The task of finding the number of equivalence classes for all integers 
n and /c, 1 < < n is an open scientific problem. We partially solve 
this problem by making a computer program to find this number for some 
(not great) values of n and k. Moreover, using bitwise operations, our 
algorithm will receive one representative from each equivalence class without 
examining the whole set A^. 
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3 Bitwise operations 



Bitwise operations can be applied for integer data type only, i.e. they cannot 
be used for float and double types. For the definition of the bitwise opera- 
tions and some of their elementary applications could be seen, for example, 
in P El El E]. 

We assume, as usual that bits numbering in variables starts from right 
to left, and that the number of the very right one is 0. 

Let x,y and z are integer variables or constants of one type, for which 
bits are needed. Let x and y are initialized (if they are variables) and let the 
assignment z = x & y; {bitwise AND), or z = x I y; {bitwise inclusive 
OR), or z = X " y; {bitwise exclusive OR), or z = ~x; {bitwise NOT) be 
made. For each i = 0,1,2, . . . ,w — 1, the new contents of the i-th bit in z 
will be as it is presented in the Table [TJ 
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Table 1: Bitwise operations 



In case that k is a nonnegative integer, then the statement z = x<<k 
{bitwise shift left) will write (i + k) in the bit of z the value of the k bit of 
X, where i = 0,1, . . . ,w — k — 1, and the very right k bits of x will be filled 
by zeroes. This operation is equivalent to a multiplication of x by 2'^. 

The statement z=x>>k {bitwise shift right) works the similar way. But 
we must be careful if we use the programming language C or C + +, as 
in various programming environments (see for example in [TU]) this opera- 
tion has different interpretations - somewhere k bits of z from the very left 
place are compulsory filled by (logical displacement), and elsewhere the 
very left k bits of z are filled with the value from the very left (sign) bit; 
i.e. if the number is negative, then the filling will be with 1 (arithmetic dis- 
placement). Therefore it is recommended to use unsigned type of variables 
(if the opposite is not necessary) while working with bitwise operations. In 
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the Java programming language, this problem is solved by introducing the 
two different operators: z=x>>k and z=x>»k p]. 
Bitwise operations are left associative. 

The priority of operations in descending order is as follows: ~ (bitwise 
NOT); the arithmetic operations * (multiply), / (divide), % (remainder 
or modulus); the arithmetic operations + (addition) - (subtraction); the 
bitwise operations << and >>; the relational operations <, >, <=, >=, 
==, !=; the bitwise operations and |; the logical operations && and ||. 



To compute the value of the z-th bit of an integer variable x we can use 
the function: 

int BitValue(int x, unsigned int i) { 
return ( (x & l«i) == ) ? : 1; 
} 

The next function prints an integer in binary notation. We don't con- 
sider and we don't print the sign of integer. For this reason we work with 
\n\ 

void DecToBin(int n) 
{ 

n = abs (n) ; 
int b; 

int d = sizeof (int) *8 - 1; 

while ( d>0 && (n & l«(d-l) ) == ) d--; 

while (d>=0) 



{ 

b= l«(d-l) & n ? 1 : 0; 

cout<<b; 

d-; 

} 



The following function calculates the number of I's in the binary repre- 
sentation of an integer n. Again we ignore the sign of the number. 

int NumbOf_l(int n) 

-c 

n = abs (n) ; 



} 
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int temp=0; 

int d = sizeof (int) *8 - 1; 
for (int i=0; i<d; i++) 

if (n & l«i) temp++; 
return temp; 

} 

4 Description and implementation of the algo- 
rithm 

Let N be the set of natural numbers and let 

Tn = {{xi,X2, . . . ,Xn) \ Xi eN, i = 1,2,. . . ,n} 
An one to one corresponding 

LP : Bn ^ Tn 

which is based on the binary presentation of the natural numbers, is de- 
scribed in [6]. li A ^ Bn and f{A) = {xi, X2, ■ ■ ■ Xn), then i-th row of A is 
integer Xi written in binary notation. 

In [3], it is proved that the representation of the elements of Bn us- 
ing ordered n-tuples of natural numbers leads to making a fast and saving 
memory algorithms. 

Let A ^ Bn and let x = {xi,X2, ■ ■ ■ , Xn) = '^{A). Then we denote 

X* = ^{A'), 

where A* G Bn is the transpose of the matrix A . 

Let X = (xi, X2, . . . , x„) and let x* = {yi,y2, . . . ,yn)- x we will call 
canonical element, ii Xi < X2 < ■ ■ ■ < Xn and yi < 22 <■■■< yn- 

Proposition 1 There is un unique canonical element in every equivalence 
class of factor-set I . 

The proof of proposition [T] is within the reach of any student who has 
successfully studied the properties of the binary system concept and we will 
miss it here. 
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Proposition 1 is the base of our algorithm, which we describe in brief 
below. For its implementation, we will use also the functions shown in 
section 3. 

As it is well known, there are exactly 2" nonnegative integers, which are 
presented with no more than n digits in binary notation. We need to select 
all of them, which have exactly k I's in binary notation. Their number 
is (^) <^ 2". We could use the function NumbOf _1 (int) from section [31 
but then we have to use it for each integer from the interval [0,2" — 1], 
i.e. 2" times. We will describe an algorithm that directly receives the 
necessary elements without checking whether any integer m G [0, 2" — 1] 
satisfies the conditions. We will remember the result in the array p [] of 
size c = (^). Moreover, the obtained array is sorted in ascending order 
and there are no duplicate elements. The algorithm is based on the fact 
that the set of all ordered m-tuples B"^ = 62, . . . , bm), hi (z B = {0, 1}, 
i = 1,2, ... ,771, m = 1,2, ... ,71, is partitioned into two disjoint subsets 
B"' = M1UM2, MinM2 = 0, where Mi = {(61,62, ■ ■ ■ , 6m) I 61 = 0} and 
A^2 = {(61, 62, ... , bm) I bi = 1}. The described recursive algorithm again 
uses bitwise operations. 

void DataNumb(int p[], unsigned int n, unsigned int k, int& c) 
{ 

if (k==0) 
{ 

c = 1; 
p[0] = 0; 
} 

else if (k==n) 
i 

c = 1; 
p[0] = 0; 

for (int i=0; i<k; i++) p[0] 1= l«i; 
} 

else 

{ 

int pi [10000] , p2 [10000] ; 
int cl, c2; 

DataNumbCpl, n-1, k, cl) ; 
DataNumb(p2, n-1, k-1, c2) ; 
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c = cl+c2; 

for (int i=0; i<cl; i++) p[i] = pl[i]; 

for (int i=0; i<c2; i++) p[cl+i] = p2 [i] I l«(n-l); 

} 



We also will use bitwise operations in constructing the next two func- 
tions. 

The function int n_tuple (int [] , int, int, int) gets all t = ("^^~^) 
(combinations with repetitions) ordered n-tuples {xi,X2, • • • , x„), where < 
Xi < X2 < ■ ■ ■ < Xn < c, Xi, i = 1,2, . . . ,n are elements of sorted array p [] 
of size c. As a result, the function returns the number of canonical elements. 

The function bool check (int [] , int) refers to the use of each received 
n-tuples. It examines whether this is a canonical element and prints it. 

bool check(int x[], int n, int k) 
{ 

int yj ; // the integer representing column (n-j) 
int yO=0; // integer preceding column j 
int b; 

for (int j=n-l; j>=0; j--) 
{ 

yj=o; 

for (int i=0; i<n; i++) 
{ 

b = l«j & x[i] ? 1 : 0; 
yj 1= b « (n-l-i); 

} 



if (yj<yO II (NumbOf.Kyj) != k)) return false; 

yo = yj; 

} 

//We have received a canonical element. Print it: 
for (int i=0; i<n; i++) cout«x[i]«" "; 
cout«'\n' ; 
return true; 
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int n_tuple(int p[], int n, int k, int c) 
{ 

int t=0; 

int a [n] , x [n] ; 

int indx = n-1; 

for (int i=0; i<n; i++) a[i]=0; 

while (indx >= 0) 

{ 

for (int i=indx+l; i<n; i++) a[i] = a [indx] ; 

for (int i=0; i<n; i++) x[i] = p[a[i]]; 
if (check(x,n,k) ) t++; 

indx = n-1; 
a [indx] ++; 

while (indx>=0 && a[indx]==c) 
{ 

indx-- ; 
a [indx] ++; 

} 

} 

return t ; 

} 

The description of the main function, we leave to the reader. 
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Table 2: The number of equivalence classes for 1 < A; < n < 8 
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5 Conclusion 



The number of equivalence classes for 1 < A; < n < 8 are given in Table 
m which is obtained through the work of the algorithms described in this 
paper. 

The ideas described in this article can be used for finding the cardinality 
of other factor-sets of binary matrices 
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